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. . - Abstract 

j^i^ . Kerr-Schild solutions to the vacuum Einstein equations are consid- 

' ered from the viewpoint of integral equations. We show that, for a class 

, of Kerr-Schild fields, the stress-energy tensor can be regarded as a total 

, ^ ■ divergence in Minkowski spacetime. If one assumes that Minkowski co- 

^ l' ordinates cover the entire manifold (no maximal extension), then Gauss' 

, theorem can be used to reveal the nature of any sources present. For the 

Schwarzschild and Vaidya solutions the fields are shown to result from a 
(5-function point source. For the Reissner-Nordstrom solution we find that 

' . inclusion of the gravitational fields removes the divergent self-energy famil- 
iar from classical electromagnetism. For more general solutions a complex 

^ ^ structure is seen to arise in a natural, geometric manner with the role of 

I . the unit imaginary fulfilled by the spacetime pseudoscalar. The Kerr solu- 

QQ ' tion is analysed leading to a novel picture of its global properties. Gauss' 

, theorem reveals the presence of a disk of tension surrounded by the mat- 

' ter ring singularity. Remarkably, the tension profile over this disk has a 

' simple classical interpretation. It is also shown that the matter in the ring 

, follows a light-like path, as one expects for the endpoint of rotating, col- 

' lapsing matter. Some implications of these results for physically-realistic 

^ ' black holes are discussed. 

7'. 

1 Introduction 

, ^ ' Many of the important solutions to the Einstein field equations can be ex- 

. pressed in Kerr-Schild form (see, for example, the discussion in ^). These 

$H ' include all black hole solutions, and a range of solutions representing radiation. 

Here we analyse solutions of Kerr-Schild type from the viewpoint of the gauge 
theory approach to gravity |31 E] . In this approach the gravitational fields 
are gauge fields defined over a flat Minkowski spacetime. These fields ensure 
that all relations between physical quantities are independent of the position 
and orientation of the matter fields — a scheme that ensures that the back- 
ground spacetime plays no dynamic role in the physics and has no measurable 
properties. Kerr-Schild metrics are constructed from a null vector field in the 
background Minkowski spacetime, so are particularly well-suited to analysis via 
this gauge-theoretic approach. In this paper we show that, for all fields of 
Kerr-Schild type, the Einstein tensor is a total divergence in the background 
Minkowski spacetime. Various consequences of this result are explored. Gauss' 
theorem is used to convert volume integrals of the Einstein tensor to surface 
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integrals, enabling us to probe the nature of the matter singularities generating 
the gravitational fields. 

The gauge-theory viewpoint always produces a metric that satisfies the Ein- 
stein equations (or their generalisation to include torsion). But working with 
fields defined over a Minkowski background does place additional restrictions on 
the form of the solutions. For example, general relativity admits two possibilities 
when dealing with the Kerr solution jS] 

1. The complete Kerr manifold can be covered by a single set of Minkowski 
coordinates. This implies a discontinuity in the fields over the entire disk 
region bounded by the matter singularity. 

2. The fields arc smooth everywhere away from the ring, but an observer 
passing through the ring emerges in a new, asymptotically flat, region. 
This is achieved by extending the radial coordinate r to negative values, 
producing the maximally-extended Kerr spacetime. 

By adopting a flat-spacetime, gauge-theory formulation we restrict ourselves 
to considering case 1 only. This can be justified on the grounds that the 
full, maximally-extended Kerr solution is not thought to be a feasible endpoint 
for any collapse process. Similar comments apply to the maximally-extended 
Schwarzschild and Reissner-Nordstrom solutions, neither of which are consid- 
ered here. 

The first applications we consider are to spherically-symmetric fields, con- 
centrating on the Schwarzschild, Reissner-Nordstrom and Vaidya solutions. In 
all cases the integrals provide sensible results for the total energy contained in 
the fields, with the mass contribution to the energy residing in a point-source 
(5-function. For the Reissner-Nordstrom solution the inclusion of gravitational 
fields removes the infinite electromagnetic self-energy for a point charge familiar 
from classical electromagnetism This result is achieved without requiring 
any form of regularisation procedure, and ensures that the total electromagnetic 
self-energy is zero. 

We next turn to more general fields following the work of Schiffer et al. [Jj. 
These authors showed that stationary Kerr-Schild vacuum solutions are gener- 
ated by a single, complex generating function. This complex structure underlies 
the 'trick' by which the Kerr solution is obtained from the Schwarzschild so- 
lution via a complex 'coordinate transformation' (This is a trick because 
there is no a priori justification for expecting the complex transformation to re- 
sult in a new vacuum solution.) The complex structure associated with vacuum 
Kerr-Schild fields is shown here to have a simple geometric origin, with the role 
of the unit imaginary fulfilled by the spacetime pseudoscalar — the same entity 
that is responsible for duality transformations of the Riemann tensor. 

The remainder of this paper deals with a detailed analysis of the Kerr solu- 
tion. For this we require a careful choice of branch cut in the complex square 
route in the generating function. Once this is made. Gauss' theorem reveals the 
detailed structure of the singular region, confirming that the matter is concen- 
trated in a ring that circulates on a lightlike trajectory. This is as one would 
expect, since the Riemann tensor only diverges on a ring, and special relativity 
alone is sufficient to predict that rotating collapsing matter will fall inwards 
until its velocity becomes lightlike. A more surprising result is obtained from 
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considering integrals inside the ring, which reveal the presence of a disc of pla- 
nar tension [5]. This tension is isotropic over the disk and has a simple radial 
dependence, rising to oo at the ring. Remarkably, the functional form of the ten- 
sion has a simple non-gravitational interpretation. In non-relativistic dynamics 
a membrane holding together a rotating ring of disconnected particles would 
be under a constant tension. When spccial-rclativistic effects are included the 
picture is altered by the fact that tension can act as a source of inertia. This 
introduces a radial dependence into the tension, the functional form of which is 
precisely that which lies at the heart of the Kerr solution. These conclusions are 
gauge invariant and are not artifacts of the use of the background spacetime. 
This is demonstrated by eigen-decompositions of the stress-energy and Riemann 
tensors, from which we extract the gauge-invariant information. 

There has been considerable debate over many years surrounding the nature 
of sources for the Kerr metric. Many physicists have attempted to construct 
extended sources for which the Kerr metric could represent the external geom- 
etry (see Krasinski for an early review ^Hl)- More recently, a series of authors 
have constructed disk sources for the Kerr metric ^2^]^]. These solutions 
represent extended sources and do not have horizons present. The present work 
is of a different nature, dealing solely with the structure of the singular region 
— the endpoint of a collapse process. The first authors to consider this were 
Newman & Janis |H| and Isreal We disagree with Isreal's result for the en- 
ergy distribution over the disk, agreeing instead with Hamity's later result 
Our techniques enable us to go some way beyond Hamity's description, both in 
revealing the physical properties of the disk and in understanding the nature of 
the singularity around the ring. The simple physics of the disk was first pointed 
out in [H]. 

Many of the calculations here are simplified by using the language of 'space- 
time algebra' [3 1151 [TCj . This is crucial to understanding the geometric nature 
of the complex structure at the heart of Kcrr-Schild solutions. The algebraic 
structure of the spacetime algebra is that of the Dirac 7-matrices. Using this 
algebraic structure one can develop a mathematical language that is adept at 
describing many aspects of relativistic physics. This language includes a calcu- 
lus that is somewhat more powerful than any available in alternative languages. 
The gauge theory of gravity developed in |2] takes on its most natural and 
compelling form when expressed in the spacetime algebra. We start with an 
introduction to the spacetime algebra, giving the necessary conventions and 
notations. Further details can be found in |21 El IIM IT7j and references con- 
tained therein. Reference |21 includes an appendix describing how to convert 
between spacetime algebra and more conventional tensor calculus. Natural units 
(G = c = Co = 1) s-rc employed throughout this paper. 

1.1 Spacetime algebra 

The basic algebraic structure behind the spacetime algebra will be familiar 
to most physicists in the guise of the algebra of the Dirac 7-matrices. The 
geometric interpretation the spacetime algebra attaches to this algebra may be 
less familiar, though it is remarkably well-suited to most problems in relativistic 
physics 13 El El- The spacetime algebra is generated by four vectors {7^}, IJ, = 
. . . 3, equipped with an associative (Clifford) product denoted by juxtaposition. 
The symmetric and antisymmetric parts of this product define the inner and 
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outer products, and are denoted with a dot and a wedge respectively, so 



lt,-lu = lili^lu + 71/ 7m) = 'Hi^u = diag(+ 



(1.1) 



and 



(1.2) 



The outer product of two vectors defines a bivector — a directed plane segment 
representing the plane defined by the two vectors. A full basis for the spacetime 
algebra is provided by the set 
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{7m} {o-kJtTk} {/7m} 



1 scalar 4 vectors 6 bivectors 4 trivectors 1 pseudoscalar, 
grade grade 1 grade 2 grade 3 grade 4 



where 



and 



k = l...i 



70 71 72 73 



cr icr2cr3. 



(1.3) 



(1.4) 



(1.5) 



The pseudoscalar / squares to —1, anticommutes with all odd-grade elements 
and commutes with even grade elements. Both the {<Tfc} and {7^} are algebraic 
entities with clear geometric significance. They should not be thought of as 
matrices acting on an internal spin space. (The same symbols as employed in 
quantum theory are used here simply because the algebraic relations are the 
same.) 

An arbitrary real superposition of the basis elements 1)1.3(1 is called a 'multi- 
vector' and these inherit the associative Clifford product of the {7^} generators. 
The inner and outer products with a vector a are of particular importance. For 
these we write 



(1.6) 



The outer and geometric products are associative, but the inner product is not. 
We also employ the commutator product, 



Ay,B^ UAB- BA). 



(1.7) 



Vectors are usually denoted in lower case Latin, a = a^7^, or Greek for basis 
frame vectors. In the absence of brackets the inner, outer and commutator 
products take precedence over geometric products. 

An inertial system is picked out by a future-pointing timelike (unit) vector. If 
this is chosen to be the 70 direction then the 70-vcctor determines a map between 
spacetime vectors a = a^^^ and the even subalgebra of the full spacetime algebra 
via 

a-fo = ao + a, (1.8) 



where 



and 



aA7o. 



(1.9) 



The 'relative vector' a can be decomposed in the {crk} frame and represents 
a spatial vector as seen by an observer in the 70-frame. Relative (or spatial) 
vectors in the 70-system are written in bold type to record the fact that they 
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are in fact spacetime bivectors. This distinguishes them from spacetime vectors, 
which are left in normal type. The {cTk} generate the (Pauli) algebra of three- 
dimensional space, and we occasionally require that the dot and wedge symbols 
define the three-dimensional inner and outer products. The convention we adopt 
is that, if both arguments of a dot or wedge product are written in bold, then 
the product takes its three-dimensional meaning. For example, aAb is a relative 
bivector, and so also a spacetime bivector, and not a spacetime four-vector. 
The vector derivative, V, is defined by 

v = r^ (1.10) 

where the {x'^} are a set of Cartesian coordinates and the {7''} are the reciprocal 
frame to the associated coordinate frame {7^}, i.e. ^^-^u = Sf^. The spacetime 
split of the vector derivative V goes through slightly differently, since we require 
that the V symbol agrees with its conventional three-dimensional meaning. This 
is achieved by writing 

7oV = Ot + V, (1.11) 

so that V = <Jidj.i. The V operator has the algebraic properties of a vector, 
and often acts on objects to which it is not adjacent. The 'overdot' notation is 
a convenient way to encode this: 

V^^ = 7^A^. (1.12) 

The V operator acts on the object to its immediate right unless brackets or 
overdots are present. If brackets are present then V operates on everything in 
the bracket, so that, for example, V{AB) ~ V AB + V AI3. The same rules 
apply to V. 

One of the two gravitational gauge fields is the (position-dependent) linear 
function h(a), which maps vectors to vectors (where a is the vector argument). 
Linear functions of this type have their action extended to general multivectors 
via the rule 

h(aA6---Ac) = h(a) Ah(6) A- • • Ah(c), (1.13) 

which defines a grade-preserving linear operation. The pseudoscalar is unique 
up to a scale factor, and the determinant is defined by 

h(/) = dct(h)/. (1.14) 

The adjoint is denoted with an overbar, h(a). The fimction h(a) and its adjoint 
are related by ^Hl 

A-h(S,.) h(h(Ar)-Bs) r<s 

h{Ar)-Bs^h{Ar-h{Bs)) r > s. (1.15) 

A number of manipulations in linear algebra are simplified by using the 
vector derivative in place of frame contractions. For example, the trace of h(a) 
can be written as 

Tr(h)=7^-h(7M) = 9a-h(a), (1.16) 
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where da is the vector derivative with respect to a. The following results are 
also useful: 

daa-Ar = rAr (1-17) 
daaAAr ^ {n- r)Ar (1.18) 
daAra = {-lY{n - 2r)Ar, (1.19) 

where Ar is a multivector of grade r and n is the dimension of the space. 
1.2 The field equations 

The gravitational gauge fields are a linear function h(a) mapping vectors to 
vectors and a linear function r2(a) mapping vectors to bivectors. Both of these 
gauge fields have an arbitrary position dependence. The gauge-theoretic origin 
of these fields is described in [SJ 0] . The gauge fields are related by the equation 

2n{a) ^ -h(VAg(a)) + h-'^{db)A{a-Vh{b)), (1.20) 

where 

g(a) = h-ih-i(a). (1.21) 

The argument of the linear function, usually denoted by a vector a or b, is 
always assumed to be independent of position. To recover the more conventional 
representation of general relativity we introduce an arbitrary set of coordinates 
a;'', with the associated coordinate frame vectors, 




(1.22) 



With e'^ denoting the reciprocal frame vectors we then define the vectors 

5^ = h-i(e^), g^ = Ken. (1.23) 

In terms of these the metric tensor is defined by 

9t,u^9t,-9v (1.24) 

The h(a) field ensures that one only ever has to make 'fiat-space' contractions, 
which is an attractive feature of the gauge-theory approach. 

The field strength corresponding to the 17 (a) gauge field is defined by 

R(aA6) = a-Vn{b) - b-VVL{a) + n{a) x VL{b) (1.25) 

and is a linear function mapping bivectors to bivectors. From this the covariant 
Riemann tensor is defined by 

7^(aA&) = R(h(aA5)). (1.26) 

We often write this in the form TZ(B), where B is an arbitrary (constant) bivcc- 
tor argument. The tensor components of the Riemann tensor arc recovered by 
writing 

R'^.pa = i9''Ag,)-nig,Agp). (1.27) 
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The Ricci and Einstein tensors are defined from tlie Riemann tensor in tlie 
obvious way, 

Ricci Tensor: 7^(6) = •7^(aA6) (1.28) 
Ricci Scalar: TZ ^ da-n{a) (1.29) 
Einstein Tensor: ^(a) = 7^(a) - ia7^. (1.30) 

Again, the tensor components of the Ricci and Einstein tensors are easily re- 
covered. 

1.3 Kerr-Schild fields 

We are interested in fields of the form 

h{a)^a + a-ll (1.31) 

where I is a (fiat-space) null vector, P = 0. This is the gauge theory analogue 
of the Kerr-Schild ansatz. The function H1.31(l extends to act on multivectors 
as 

h{A) = h{A) ^ A + A-ll, (1.32) 
and we see immediately that dct(h) = 1. The following results are also useful: 

h-\A) ^h-\A) ^ A- A-ll (1.33) 
g{A) = A-2A-ll (1.34) 
h{l)^l. (1.35) 

In terms of an orthonormal coordinate frame 7^ we can write 

gti^lii-lpil (1-36) 

which confirms that the metric is given by 

9tiv = Vnu - '^Ifj.lu, (1-37) 

where Ty^j/ is the flat Minkowski metric tensor. 

The fl{a) field defined by (|1.31|1 has the simple form 

n{a) = h(VA(a-n)) 

= \/A{a-ll)-a-lvAl (1.38) 

where 

v^l-Vl. (1.39) 
It follows from the fact that / is null that 

l-v = Q (1.40) 

and 

17(0=0. (1.41) 
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Following the route adopted by Chandrasekhar ^1 Section 57], we next 
form the quantity 

/•7^(0 = l-{da-n{aM)) 
= {lAda)-R{aAl) 

= {lAda)-{a-<7n{l) l-Wn{a)). (1.42) 

Substituting equation 1)1.38(1 into the above we find that 

l-TZ{l) ^ {lAda)-{-<7{{a-Wl)-i)M-l-WWA{a-ll)) 
= da-l{a-Vl)-v~l-V(y-ll + v)-l 
= v^-{hVv)-l 

= 2v^. (1.43) 

If we are looking solely for vacuum solutions, then we can conclude from this 
that V must be null. Since = 0, it follows that v must be parallel to /, 

V = (jyl, (1.44) 

where is a scalar field. We will restrict attention to solutions for which this re- 
lation does hold, even if matter is present. (This places a restriction on the form 
of matter distributions that we can consider.) It follows from equation ((1.44() 
that ri(a) reduces to the simpler form 

n{a) h{a-ll). (1.45) 

The Riemann tensor now splits into terms that are second-order and fourth- 
order in I. The fourth-order contribution is 

7^4(aA6) = -tl{{{ahb)-ll)-V) + VL{a)xVL{h). (1.46) 

After some rearrangement this can be brought to the form 

7^4(B) = \da-db{a-Wl)lBl{h-^l) - \{a-Vl)-{h-^l) dJBldb. (1.47) 

Both the contraction, da ■ 'R-{a A 6), and the protraction, da A Tl{a A fe), of this 
contribution to the Riemann tensor vanish. This can be seen from the result 
that 

daFiaAbF2 = daFiiab - a-b)F2 = -feFiFa, (1.48) 

which holds for any two bivectors Fi and F2. The presence of the null vector I 
in the analogous terms in TZi{B) ensures that 

dani{aAb)=0, (1.49) 

so that TZ4{B) makes no contribution to the Ricci tensor. 

The only part of TZ{B) that contributes to the Einstein tensor is therefore 
the second-order term 

n2{aAb) = a-Vn{b) ~b-Vn{a). (1.50) 
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Contracting this and setting the result to zero we find that the vacuum Einstein 
equations reduce to solving the equation 

n{a) = V -Qia) - a-V db-n{b) = 0. (1.51) 

The Ricci scalar and Einstein tensor are now straightforward to calculate: 

n = -2V-{da-n{a)) (1.52) 

and 

Gia) = V-{n{a) - aA{db-n{b))). (1.53) 

The formulae for ri(a) H1.45|l and G{a) are valid for any Kerr-Schild type so- 
lution for which ^ • = (pi. For such fields the Einstein tensor l|1.53|l is a 
total divergence in Minkowski spacetime. In general, the field equations will be 
satisfied everywhere except for some singular region over which the fields are 
discontinuous. This singular region contains the source of the fields. In this 
paper we assume that the entire solution to the Einstein equations is described 
by fields defined over a single Minkowski spacetime, so that the manifold has not 
been subjected to maximal extension. In this case we can use Gauss' theorem 
straightforwardly to convert volume integrals over the source region to surface 
integrals and so learn how the source matter is distributed. For the case of static 
fields, Virbhadra gave a formula which agrees with H1.53|l for the timelike 
component a = 70, but the fact that the expression is a total divergence was 
not exploited. 



2 Spherically-symmetric solutions 

As our first application we consider spherically-symmetric solutions. For these 
it is useful to introduce a set of polar coordinates: 



t = x-^o cos9 = x-7 /r 



(2.1) 



T = \J (a;A7o)^ tanc/) = (a;-7^)/(x-7"'"). 

We also define 

= a::A7o7o/7', cr,. = 6^70, (2.2) 

and 

e± = 70 ± e^. (2.3) 
For spherically-symmetric solutions I can be written in the form 

/ = \/^e±, (2.4) 

where a' — a'{t, r). For fields of this type it is a simple matter to demonstrate 
that the fourth-order contribution to the Riemann tensor H1.47|l vanishes. To 
see this consider the case of e_)_, for which we obtain 

/2 

7^4(B) = ^{-da-dba-V(Tr{l - (Tr)B{l + (Tr)b-V(Tr 

+ {a-\/crr)-{b-\7crr) daC+Be+db) 

= — (V(l - Cr,.)B(l + (Tr)&r - V(l - Cr,,)B(l + (Tr)crr) 

4r ' 
= 0, (2.5) 
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with the same result holding for e_. It follows that the Riemann tensor is given 
entirely by (|1.50|l , which is also a total divergence and so can be analysed using 
Gauss' theorem. We now turn to three applications of these results. 



2.1 The Schwarzschild solution 

The simplest solution to the field equations is the Schwarzschild solution, ob- 
tained from 

a' = M/r, ? = (70 - e,.), (2.6) 

as we confirm in section 13 The line element generated by this solution is that 
of the advanced Eddington-Finkelstein form of the Schwarzschild solution. The 
Riemann tensor for the solution 1)2. 6|l can be constructed using equation H1.50|l . 
from which we find 

7^(a) = a-Vfi(7o) 

= a.V(VA(M(7o-e,)/r)) 

= Ma-Vj, (2.7) 



and 



7^(/6) = 0((/6).V7o) 



VA /6A(TrAV<Tr7o 
^^^^.(bAj). (2.8) 



Away from the origin, these derivatives evaluate to 



and 



TZ{a) = ^{a — ia-(Tr (Tr), (2.9) 



IM 

n{Ib) ^ —{b-3b-(Tr(Tr), (2.10) 



so we can write the vacuum Riemann tensor in the form 



M 

n{B) = - — {B + 2,(TrBcTr). (2.11) 



Self duality of the vacuum Riemann tensor has the simple expression TZ(IB) — 
nZ{B) in the spacetime algebra formalism [2]). The form of equation (|2.11|l 
shows that the Riemann tensor is manifestly self-dual. This form of the Riemann 
tensor (|2.11() for the Schwarzschild solution was first given in [50] and . 

The form of the Riemann tensor in equation H2.11|l is valid everywhere away 
from the singularity. To study the form of the singularity, we return to the 
diiferential expressions for the Riemann tensor and integrate over a sphere of 
radius tq, centered on the origin. Using Gauss' theorem to convert the volume 
integrals to surface integrals, we obtain 

/ (fxn{a)^M d(l) ^6* sine a- cr^cTr = ——a, (2.12) 

Jr<ra Jo Jo 3 
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and 



\ d^xTZilb)^ d(p d9 sine I(Tr-{bA(Tr)^-^^^^ lb. (2.13) 

Jr<ra Jo Jo 3 

These results combine to give 

[ d^xn{B) = 1^(5-7070 - 2BA7070) 

Jr<ro 

- -^(S + 37oB7o), (2.14) 

which contracts to yield 

d^xn{a) = 47rM 70070 (2.15) 

d^xn = -8TTM (2.16) 

d^xGia) = 87rMa-7o7o. (2.17) 



Since TZ{a) — everywhere except for the origin, the integrals ifTT^ - lfnTjl can 
be taken over any region of space enclosing the origin. It is clear then that 
the solution represents a point source of matter, with the matter stress-energy 
tensor given by 

T(a) = A/,5(a;)a-7o7o- (2.18) 

The same conclusion was reached in j2|, where the calculations were performed 
in a different gauge. This result confirms Feynman's speculation in Lecture 15 
of 1221 that "it will not be possible to demonstrate that G^^ = everywhere, 
but rather that G^,y = S{x), or something of the kind" . The integrals performed 
above are not gauge invariant, but gauge- invariant information is extracted from 
them in the form of the matter stress-energy tensor H2.18|l . Furthermore, the 
integral of the Ricci scalar provides a direct measure of the mass of the source, 
without the need to resort to constructing integrals in an asymptotically flat 
region of spacetime. 

2.2 The Reissner-Nordstrom solution 

The Reissner-Nordstrom solution can be written in the form 

h(a) = a + 77a-e_ e_ (2-19) 
where q the charge of the source (in natural units) and 

The Einstein tensor for this solution is 

g{a) = V-(V(r/a-e_)Ae_ - V- (?7e_) aAe_) . (2.21) 
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Away from the origin we know that the mass term can be ignored, which leaves 



S(7.) = |^V.(v(f:i;i:i)Ae_) = -^..rr-,,^.. (2.23) 

These combine to give a corresponding matter stress-energy tensor of 

^(«) = ^^(a) = -^^a^ (2.24) 

where J- = q(Tr/{4:Trr'^). This is the expected form for the electromagnetic 
stress-energy tensor due to a point source of charge q. (See for a detailed 
explanation of how to handle electromagnetism in gauge-theory gravity.) 

To study the behaviour of the fields near the origin we return to the differ- 
ential form for G{a) and again construct integrals over a sphere of radius vq. 
For this case we find that 

d^xGi-jo)^ [ (Tr - cr^^ 70 

SttAI - ^) 70, (2.25) 



and 



/ d^xg{a-in)= I d^a;|-7ovY4raAcr,, 



= 1^070, (2.26) 
3ro 

which combine to give 

/ d3xr(a) =Ma-7o7o + :^r— ("-4a-7o7o)- (2.27) 

The mass term here is precisely as expected and shows again that a point source 
is located at the origin. The electromagnetic contribution is traceless, as one 
expects for the electromagnetic stress-energy tensor. Focusing attention on the 
7o-frame energy component of the stress-energy tensor, we see that 

d3x7o-r(7o)-M--^. (2.28) 

This result was also obtained by Tod |22|, who calculated the quasi-local mass 
for the Reissner-Nordstrom solution as defined by Penrose . Tod argued that 
this result implies that a source for the Reissner-Nordstrom solution should have 
r > q'^/{8TrM) at the surface in order to meet the dominant energy condition. 
However, this misses the point that the negative contribution to the integral 
comes entirely from the origin. Everywhere off the origin the stress-energy 
tensor satisfies the dominant energy condition. Taking the integrals over the 
volume defined by ro < r < oo we find that the electromagnetic field energy is 
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j (87rro), which agrees with the formula given by Virbhadra |19| and is simply 
the classical result. 

The electromagnetic contribution to (|2.28() is negative and finite for all fi- 
nite ro, and tends to zero as the integral extends over all space. This is in 
stark contrast to the standard picture from classical electromagnetism, where 
the integral of the 70-frame energy jl diverges for the interior of any surface 
enclosing the origin — the classical self-energy problem discussed by many au- 
thors (see |H| , for example) . Inclusion of the gravitational field has removed 
this divergence, ensuring that the total electromagnetic self-energy is zero. The 
manner in which this regularisation is achieved is both simple and instructive. 
The electromagnetic energy density is rewritten as 



2 32..^ 32.^- bJ' (2-2^) 

so that the integral over space of the electromagnetic energy density can be 
converted to a surface integral, recovering the contribution to H2.28() . Since the 
electromagnetic energy density near a point source is very large, it is unsurpris- 
ing that the inclusion of gravity has significant consequences, and these clearly 
have implications for the status of self-energies in classical field theory. However, 
since only classical fields are employed above, it is not clear whether this result 
has similar implications for the divergent self-energies encountered in QED. 

2.3 The Vaidya solution 

As a final example of the use of integral equations for spherically-symmetric 
Kerr-Schild fields, we consider Vaidya's 'shining star' solution This is gen- 
erated by the field 

h(a) =Q+ a-e+e+, (2.30) 

r 

which is clearly similar to the Schwarzschild solution, except that now the mass 
/X = [i(t — r) is variable and the null geodesies e-|_ are outgoing rather than 
incoming. The solution (|2.30|) is clearly of Kerr-Schild type, and defining I by 

; = V^e+, (2.31) 

we find that 

'■v-(^)"^.-v((^)"%.)^4(^)"'. 

so that equation H1.44|l is satisfied. The Einstein tensor for ()2.30|1 is 

g{a) = V-{^a-e^cT^ (2.33) 
and away from the origin (where we can set V-((Tr/r^) — 0) this becomes 

e(a) = -^a-e+e+, (2.34) 

where /i = 9t/^. This tensor represents a radially-symmetric flux of outgoing 
massless particles. Again, the presence of a (5-function point source at the origin 
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can be inferred from the differential form of the Einstein tensor. By evaluating 
the integral of Q{a) over a sphere centred on the origin, and shrinking the radius 
to zero, we find that 

2/i 

G{a) = --^a-e+ e+ + 87r/x5(a;)a-7o 70- (2.35) 

The solution (|2.3()|) therefore describes a point mass at rest at the origin which 
is losing mass at some arbitrary rate. This is borne out by the Riemann tensor, 

7^(B) = -^B-e+ e+ - ^(S + 3cr,,Bcr^), (2.36) 



which exhibits a neat split into a source term describing the energy outflow and 
a Weyl term due to the point mass at the origin. 

The fact that the Einstein tensor is given by the divergence of a bivector 
implies that 

V-e(a)=0. (2.37) 
We can therefore define a conserved total energy E by 

= 87r^(-oo). (2.38) 

The total conserved energy is therefore determined by the mass of the source at 
t = —00, before it began radiating, which is clearly a sensible result. A conserved 
energy of this form will exist for any Kerr-Schild field of the type H1.31|l . provided 
that the null vector I satisfies l-Vl = 61. 



3 Stationary vacuum solutions 

We now turn to a more general analysis, dropping the requirement of spherical 
symmetry. As we are ultimately interested in the Kerr solution, however, we do 
restrict to stationary, vacuum solutions. For these we write I in the form 

l^V^n (3.1) 

where 

n = 70 - n7o, (3.2) 

= 1, and a' and n arc functions of the spatial position vector x = xAjo only. 
This is the most general form for a stationary, Kerr-Schild field. The condition 
that Z-VZ (pi immediately yields 

-n- V(Va'(7o - n-fo)) = ^(70 - ^70) (3.3) 

hence 

(j>=-n-VV^, (3.4) 

and 

n-Vn = 0. (3.5) 
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The final equation shows that the integral curves of n are straight lines. These 
define possible incoming photon trajectories in space. The fact that these lines 
a straight in the background space is a gauge-specific statement, and does not 
correspond to a physically-observable property. 

For stationary vacuum fields the equation TZ(a) ~ splits into the pair of 
equations 

V-17(7o) = (3.6) 

and 

V-rj(7.,) +cr,-V(V-(a'n)n) =0. (3.7) 
To simplify equation 1)3.6(1 we need the result that 

17(7o) = VA(a'n) = -Va' - VA(a'n), (3.8) 

On splitting into spatial vector and bivector parts equation ((3.6|l reduces to 

V^a' = (3.9) 

and 

V-(VA(a'n)) = 0. (3.10) 
3.1 The hidden complex structure 

The content of the second field equation (|3.7() is summarised neatly in the equa- 
tion 

2a' 

a- Vn + V(a-n) = — — -(V(a-n))- Vn. (3.11) 

V • (a'n) 

In [Jj the authors showed that this equation implies that we can write 

a'Vn = aa/\nn — Ij3a/\n, (3.12) 
where a and /3 are two new real scalar functions. We see immediately that 

V-n = 2Q, VAn = -2//3n, (3.13) 

and it follows that 

V-(/3n) = 0, =^ n-V/3 = -2a/3. (3.14) 

Now, setting a = V in equation (|3.12|) . we obtain 

V^n = Va-V-{an)n- IVA{f3n). (3.15) 

Similarly, writing equations ((3.13(1 in the form Vn = 2(a ~ I(3ri) and differen- 
tiating we obtain 

V^n = 2Va-2IV{l3n). (3.16) 
On combining these equations we find that 

r^•Va = /32-a^ (3.17) 

and we can therefore write 

Va - IV (in = -{a^ + (P)n + 21(5{a - I(3n). (3.18) 
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This equation is simplified by employing the idempotent element N, 



iV= in7o = i(l-n), (3.19) 

which satisfies 

= N = -nN = -Nn. (3.20) 

(An idempotent is a mixed-grade multivector that squares to give itself. Space- 
time idempotents are usually closely related to null vectors.) On postmultiplying 
by N, equation l|3.18(l yields 

V7iV ^ j'^N, (3.21) 

where 

-f = a + I[3. (3.22) 
It follows that (V7)^ — 7^ so that, if we define lo by 

oj = -, (3.23) 

7 

then u) must satisfy 

[Vuf = 1. (3.24) 

This is the first of the pair of complex equations found in [7] . The novel feature 
of the derivation presented here is that the 'complex' quantity 7 is of the form 
of a scalar -I- pseudoscalar. This gives a clear geometric origin to the complex 
structure at the heart of the Kerr solution. This complex structure carries 
through to the form of the Riemann tensor, and hence to all of the observable 
quantities associated with the solution. 

Use of the idempotent N simplifies many expressions and derivations. For 
example, differentiating l|3.21|l and pre- and post-multiplying by A'^ yields 

V^7iV = ~ N\{-i^Vn-V{V-f)ri)N (3.25) 

But we know that n-V7 = —7^ and ^nN = so we can rearrange the 

final term as follows: 

NV{Vj)nN = -NnV(yj)nN 
= NV{n'V^)nN 

= -2N-f'^VnN - NVVjhN 

= -2j^N. (3.26) 

This type of rearrangement is typical of the way that one can take advantage 
of the properties of idempotents in the spacetime algebra. On substituting this 
result into equation H3.25|l we now find that 

V^jN = 0, (3.27) 

and hence 

VS = 0. (3.28) 

This is the second of the pair of complex equations found in [7]. Solving the field 
equations now reduces to finding a complex harmonic function 7 whose inverse 
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uj satisfies (Vtj)^ = 1. The above derivation reveals the geometric origin of 
this complex structure, as well as demonstrating the role of the null vector n 
through the idempotent N. 

To complete the solution we need to find forms for a' and n. For the former 
we note that 

V-(«») _ _ V-ia'n) 

2a " 2a ^ 2a' ^^^^^ 

and 

V-(VA(an)) = V-(VA(a'n)) = 0. (3.30) 

From these it is a simple matter to show that a' = Ma, where M is some 
arbitrary constant. To recover n we use equation (|3.21() in the form 

-Vw(l-n) = (1-n), Vt^*(l +n) = (1 + n) (3.31) 

where the * denotes the complex conjugation operation (which can be written 
as Lj* ~ 7oti^7o). On rearranging we obtain 

{Vuj + Vuj*)n = 2 + Vuj ~Vuj* (3.32) 

=^ = ; ,^ ^ — ^ , (3.33) 

where the angle brackets ( ) denote the projection onto the scalar part of a 
multivector. Recall here that the x symbol represents half the commutator of 
the terms on either side, and not the vector cross product. 

Some further insight into the nature of the solution and the role of the 
complex structure is obtained from the form of ri(7o). From equations (|3.13() 
to H3.16|l it is straightforward to show that 

VA(an)=/V/3. (3.34) 

It follows that f2(7o) is now given by 

17(7o) = -M{Va + VA(an)) = -MV7. (3.35) 

This shows how the harmonic function 7 generalises the scalar Newtonian po- 
tential. This gives rise to many of the novel properties of the Kerr solution. A 
further result that is useful in later calculations is that 

da-^{a) = M{a^ + 0^)n. (3.36) 
3.2 The Riemann tensor 

The Riemann tensor would be expected to contain terms of order M and A/^, but 
it is not hard to see that the latter contribution vanishes. Using equations (|1.47|l 
and H3.12|l this term can be written in the form 

A/2 

'Ri{B) = —Ia^[3{l-n)VBri{l + n). (3.37) 

But when B is the spatial bivector a we see that 

(1 - n)Vah{l + n) = (1 - n){2a-Vn - aVn)(l + n) 
= (1 - n)(27*aAn - 27*a)(l + n) 
= -27*a-n(l -n)(l + n) 
= 0, (3.38) 
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and the same result holds for B = lb. The annihilation of orthogonal idempo- 
tents in this derivation, (1 + n)(l — n) = 0, is merely a reexpression of the fact 
that n is a null vector. 

The only contribution to Tl{B) therefore comes from TZ2{B). For a spatial 
bivector this contribution can be written as 

7^(a) = 7^(aA7o) = a- V0(7o) 

= -Ma- VV7 

= -iMVaV7. (3.39) 

For a vacuum solution the Riemann tensor only contains a Weyl term, and so 
satisfies the self-duality property ^ TZ{IB) = nZ{B). It follows that we can 
write 

7^(B) = -^MVBVj, (3.40) 

for all B. This expression captures all of the terms of the Riemann tensor in 
a single, highly compact expression. Verifying that we have a vacuum solution 
now reduces to the identity 

daV ahb V7 = -bV'^-f = 0. (3.41) 



4 The Kerr solution 

The simplest solution to the pair of equations V^7 = and (Vw)^ = 1 is 
7 = 1/r. This recovers the Schwarzschild solution. To conirm this we first see 
that equation H3.33|l gives 

n^cr^, (4.1) 

which is the only possible vector consistent with spherical symmetry. The null 
vector n is given by 70 — e^, recovering the form of solution analysed in sec- 
tion l2.1l The Riemann tensor can be found directly from equation H3.4()|l . which 
yields 

-RiB) = ^MV{Bx/r^) 
M 

= -^{B + icTrBcTr) (4.2) 

recovering equation ()2.11(l . 

Since the equations V^7 ~ and (Voj)^ = 1 are invariant under complex 
'coordinate transformations', a new solution is obtained from the Schwarzschild 
solution by setting 

uj={x^ + + {z-ILfY''\ (4.3) 

This is the most general complex translation that can be applied to the Schwarz- 
schild solution 13 and generates the Kerr solution. The symbol L (L > 0) for the 
angular momentum is preferred here to the more common symbol a as we have 
already made extensive use of a as a vector variable. As was first shown in [7], 
this complex transformation justifies the 'trick' first discovered by Newman & 
Janis IHI ■ Precisely how the complex square root in (|4.3|l is defined is discussed 
further in Section ^21 



18 



From equation H3.4U|I . we can immediately construct the Riemann tensor as 
follows: 



2w3 

The spacetime bivector 

X — L Icr^ 



M fj^^.^x-LIa,^x-LIas\ ^^^^^ 



= Vuj (4.5) 



satisfies cr^^ ~ 1, so the Riemann tensor 1)4.4(1 has the same algebraic structure 
as for the Schwarzschild solution (it is type D). The only difference is that the 
eigenvalues are now complex, rather than real. The Riemann tensor is only 
singular when uj = 0, which is over the ring r = L, z = 0. This is the reason 
why the solution is conventionally referred to as containing a ring singularity. 

The structure of the fields away from the region enclosed by the ring is most 
easily seen in an oblate spheroidal coordinate system. Such a system is defined 
by: 

Lcoshit cosu = {x^ + y^Y^'^ = p (4-6) 
i sinhw sinu = z, (4-7) 

with < u < oo, — 7r/2 < v < 7r/2. These relations are summarised neatly in 
the single identity 

L cosh(M + /u) = p + Iz. (4.8) 

The basic identities for oblate spheroidal coordinates, and their relationship to 
cylindrical polar coordinates, are summarised in Tabled 

The point of adopting an oblate spheroidal coordinate system is apparent 
from the form of lo: 

w — L(cosh^w cos^u + sinh^M sin^iJ — 1 — 2/sinhu sinf)"'^^^ 

= L(sinhu — /sinw). (4-9) 

This definition of the square root ensures that u; i— > r at large distances. Equa- 
tion H3.33|l yields a unit vector n of 

2L coshu e" — ^^(coshu e" + I cosv e^)x (coshu e" — I cosv e"") 
1 + ( (coshit e" + / cosv ) (coshu e" — / cosf; e'' ) ) 

^ —^——{eu- Lcosva-A,). (4-10) 
L coshw 

As a check, 



— du 9~^<* I ftanhu cosu e„ + sinw 0-3 — (Ta\ 

L coshit Lcosh u / v coshu / 

= (4.11) 
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Cylindrical Polar Coordinates {p, (j>, z} 



p=ix^+y^y/^ 
4> ~ ta.'or^ {y I x) 

Bp = COS0 (Ti + sin(/) (T2 

60 = p(T0 ~ p{— sin^CTi + COS(/)(T2) 

epB^cTs = pi 



Oblate Spheroidal Coordinates {u, </>, v} 

L coshii cosu ~ p 
L sinhu sinu — z 

e„ = L(sinhw cosw Bp + coshu sinw ar^) 
By = L{— cosliit sinu Bp + sinhu cosu cts) 
b\ ~ b1, = L^(cosli^w — cos^w) 
e„606„ = pL^(cosli^M — cos^v)! 



Further Relations 

Bp = _L(sinhM costi 6" — coshw sinw 6") 
0-3 = L(coshw sinu 6" + sinliu cosi; 6") 
X = _L-^(sinhu coshu 6" — sinf cosu e") 



Table 1: Some basic relations for oblate spheroidal coordinates. 

and 

= 2~(cosh^u - cos^ v + cos^ v) 1, (4-12) 

cosh u 

both as required. 

The vector field n satisfies n-^n — 0, so its integral curves in flat space are 
straight lines. These can be parameterised by 

a;(A) = Lcoswq ep(0o) - A(cosi;o (T0(0o) - sin^o 0-3), (4.13) 

where L cosuq and (po are the polar coordinates for the starting point of the 
integral curve over the central disk. Plots of these integral curves are shown in 
Figure n As commented on earlier, the fact that the trajectories are straight 
lines in the background space is a feature of our chosen gauge. The same picture 
is not produced in alternative gauges, though the fact that the integral curves 
emerge from the central disk region is gauge invariant. 

4.1 Exterior integrals 

Oblate spheroidal coordinates are very useful for performing surface integrals in 
the Kerr solution over regions entirely surrounding the disk. The most conve- 
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Figure 1: Two views of the integral curves of n. The left-hand figure shows the 
view from above of a set of incoming geodesies that terminate along a diameter 
of the central disk. This pattern is rotated around the z-axis to give the full set 
of geodesies. The right-hand figure shows incoming geodesies from above and 
below the disk. These are the mirror image of each other. The ring is shown 
for clarity — it is not an integral curve of n. 



nient surfaces to consider are ellipsoids of constant u, for which the divergence 
theorem can be given in the form 



/ (fx{AVB)= I d(j) j dvpAeuB, (4.14) 

Ju'<u Jo J-Ti/2 



where A and B are general multivectors. The ^ on V indicates that the vector 
derivative acts both to the left and right, 

AVB = AVB + A{VB), (4.15) 

and the measure <fix is taken as running over u' rather than u. This slightly loose 
notation ensures that the result of the integral is a function of u. Equation (|4.14() 
accounts for all of the cases that we will encounter. 

Our aim is to explore the nature of the matter singularity through the use 
of integral equations. As a first step, we look at the total mass-energy in the 
source. Taking the surface as one of constant u we find that 

d^xg{-io)^M d(f> dwp(7oe„)-(-V7+(a' + /3')n) 

u'<u Jo J-tt/2 

Jo J~tt/2 V dv LcoshuJ 

/"^^ , /cosh^w cosw(sinh^u — sin^u) 
dv 7—T9 . 2 AO ^ 
-Tr/2 \ (smh-u + sm ij)^ 

= 87rAf7o. (4.16) 



21 



So, as the with Schwarzschild case, the total mass-energy in the 70-frame is M . 
For the spatial part we find that 

/• pliT />7r/2 

/ d^xQ{a^Q) — M 1 dcf) I du/9(7oeii)- (— V(aa-n) — VA(Q;a-nn) 

Ju'<u Jo J-tt/2 

+ (a'^ +(3^){a + aAn)). (4.17) 

Once the angular integral is performed, the 70 contribution to this integral 
becomes 

.7r/2 

27rAf7o / dv p(-du{asmv) + (a^ + coshu sinu) = 0, (4.18) 

J-7r/2 

which vanishes as the integrand is odd in v. This is reassuring, as we expect the 
integrated stress-energy tensor to be symmetric if there is no source of torsion 
hidden in the singularity. The remaining integral to be performed transforms 
to 

2tt i'Tr/2 



/ d^a;t/(a7o)7o = A/ d(j> dv pi e.^hV {(3a-n) 

Ju'<u Jo J--IT/2 



t/2 

p27T pTr/2 

dcf) j dvl{p(T^a-{dy{Pn)) - e,ja-{d^{Pn))) 

Jo J-Tr/2 
p27r pTT/2 

= M dcj) dv ip a-n{—dyp(T^ + d^By) 

Jo J-Tr/2 

= 0. (4.19) 

In terms of the corresponding matter stress-energy tensor T(a) the above results 
are summarised by 

' d^xT{a) = Ma-7o7o, (4.20) 

where the integral is over any region of space entirely enclosing the central disk. 

Integrating the matter stress-energy tensor over the entire disk region aver- 
ages out any possible angular momentum contribution. To recover the angular 
momentum we look at 

xAg{a)^t-i^Ag{a) + (a;7o)A ((70 V)- (r!(a) - aA(5b-17(5)))) . (4.21) 

The first term on the right-hand side will give zero when integrated over a region 
enclosing the disk. To simplify the remaining term we first write 

F{a)^n{a)-ah{db-VL{b)). (4.22) 

We next use the rearrangement 

(a;7o)A((7oV)-^^(a)) = (a;7o)A((7o V)-F(a))- (F(a)-7o7o + 2F(a)A7o7o), (4.23) 
to write the volume integral as 

d^xxAg{a)= d(p dvp(a;7o)A((7oe„)-F(a)) 

u'<u Jo J~it/2 

d^x (i^(a) -7070 + 2^^(a) A7o7o). (4.24) 

u' <U 
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The final volume integral involves 

n{a) - aA{db-^{b)) = MV A{aa-nn) - MaA{V ■{an)n + an-\/n) (4.25) 

which is also a total divergence and can be converted to a surface integral. For 
the 7o term we now find that 

(fxxAg{-fQ) = M 

u' <u 







/ dcp 


[ dvp( 


Jo 


J~tt/2 \ 



(a;7o)A((7oe„)-(-V7 + (a^ + I3^)n)) . (4.26) 



This simplifies down to 

f.27r /•tt/2 



M d(f> dvp{-Ix-{euAVj3) + 2aeuAn) 

Jo J-Tr/2 



/^^^ , / cosh^u cos^u(sinh^u — sin^u) 2 sinh^u cos"^w 
dv 2 ^ 2 2~ 

.Tr/2 \ (sinh u + sin u)^ sinh u + sin v, 

= -8nMLIcr3. (4.27) 
For the spatial terms we obtain 



/•2tt 




/ d4> 


[ dvpi 


Jo 


J-tt/2 \ 



d^xxAQiajo) = 

U<.Uq 

+ a;e„-(-V(aa-n) + {a^ + 0^)a) + Q;(2(e„ Aa) -nn + (e„Aa)-n) ) , (4.28) 
which reduces down to 

Ml d(f> i dv [ —2LIcr3COSva-n — 2paeu-naAn 



J-7r/2 



(4.29) 



+ pa{euAa)-n + px e„ • (— V(Q;a-n) + (a^ + /?^)a)^ 
=7\/ / (i0 / dv pi a{euAa)-n 

Jo J~it/2 \ 

+ xeu-{-V{aa-n) + [a^ + 0^)a)^. 
The final integral is best performed term by term, yielding 

J d^xxAGiji) ^ 4:nML(T2 (4.30) 

J d^xxAg{j2) =-4:nMLai (4.31) 

J d^xxAg{-/3) = 0. (4.32) 
These results combine to give 

d^xxAT{a) = ML{-a-jnI(T3 + ^ (a A70) x /erg) , (4.33) 
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where the integral is taken over any region entirely enclosing the central disk. 
This expression clearly identifies ML as the total angular momentum in the 
fields, as expected from the long-range behaviour. The expression also has 
the correct algebraic form for a symmetric stress-energy tensor. A symmetric 
stress-energy tensor has 

daST{a) = 0. (4.34) 
If this relation holds then we expect that 

<9aA(-a-7o/cr3 + i(aA7o)x(/cr3)) =0, (4.35) 

which is easily confirmed. We therefore expect that there are no sources of 
torsion hidden in the singular region. A discussion of the gauge-invariance of 
the mass-energy and angular-momentum integrals will be delayed until after we 
have a more complete understanding of the stress-energy tensor. 

4.2 The singularity 

In order to fully understand the nature of the source matter for the Kerr solution 
we must look at the region p < L. One has to be careful with the application of 
oblate spheroidal coordinates in this region, and it is safer to return to cylindrical 
polar coordinates for most calculations. Central to an understanding of this 
region is the definition of the complex square root in H4.3|) . If we consider some 
fixed p > L, then the complex function uj'^ has a real part > for all values of z 
and the square root can be defined as a smooth continuous function also with a 
real part > (see Figure|21a). This is the definition of the square root implicitly 
adopted in equation (|4.9|l with the introduction of oblate spheroidal coordinates. 
If we now consider a region where p < L and z is finite, continuity of ui requires 
that the square root still be defined to have a positive real part. This means that 
positive and negative z now correspond to different branches of the square root 
(see Figure El 6). As a result, u; is discontinuous across the entire disk p < L, 
2 = 0. This discontinuity is also easily seen in oblate spheroidal coordinates, 
for which z = 0, p < L implies that u — and sinw is discontinuous over the 
disk. The alternative, which is not considered here, is to extend the manifold so 
that passing through the disk connects an observer to a new spacetime (a new 
Riemann sheet). 

4.3 The Ricci scalar 

The simplest gauge-invariant quantity to study over the disk is the Ricci scalar, 
which is given by the total divergence 

7^ = -2V-{da-n{a)) = ~2MV-{{a^ + P^)n). (4.36) 

We compute the integral of this over an infinite cylinder centred on the z-axis 
of radius p, p < L. In converting this to a surface integral the contributions 
from the top and bottom of the cylinder can be ignored, leaving 

p /'27r poo 

/ (fxTZ=-2M d(f> dzp{a^ + P^)ep-n. (4.37) 



24 





Figure 2: The complex function ui for fixed p as a function of z. In both cases 
L = 1. The top figure is for p = 2 > L and the lower figure for p = 0.5 < L. 
The solid lines are for cu'^ and the broken lines for the square root. Continuity 
of uj for finite z requires that cu be discontiuous over the central disk. 

(As earlier, the dummy radial variable in the measure is taken as p' , so that the 
result is a function of p.) We therefore define 



/• -I /'27r /"OO 

W{p)= Sxda-T{a)/M = — d(j> dzp-f-f* Bp-n. 
Jp'<p 47r Jo J-oo 

sinhu cosu pL sinhu 



Now 



coshu L^cosh^u' 



(4.38) 
(4.39) 



and we can write 



L sinhu = M((p2 + {z- ILff^) = M . (4.40) 



We therefore only require an explicit expression for L/^ cosh u in terms of p and 
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Figure 3: Branch cuts and contours for the integral (|4.45() . The branch cuts for 
complex z shown in the top two figures follow from the global definition of lo. 
The bottom figure shows the six integrals that have to be calculated after the 
integration contour has been rotated. 



z. Such an expression is found from 

+ = L^(cosh^wcos^v + sinh^wsin^w + 1) 
= L^(cosh^w + cos^u) 



and 
so that 



^^(cosh^w - cos^O = {{p^ + - L^f + 42.2 



2L^ COSh^W = p2 ^ z2 ^ L2 + ((p2 + ^2 _ ^2)2 ^ 4^2^2^1/2^ 

On substituting these results into (|4.38(l we obtain 

R((p2 + (2 + /L)2)-l/2) 



W{p) = dzp 



(4.41) 
(4.42) 

(4.43) 
(4.44) 



p2 + Z2 + L2 + ((p2 + ^2 _ ^2)2 _^ 4^2^2)1/2 ' 

where the integrand contains a finite jump aX z — Q, {p < L) and no singularities. 
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The integral is simplified by rescaling to give 

r M((l + (z + /A)2)-i/2) 

^^'^ = y_ z^ + l + A^ + ((.^ + l-A2)2+4A^,2)i/2 ^ (4-45) 

where A ~ L/p > 1. The branch cuts for the complex square roots in this 
integral follow from the global definition of ui and arc shown in Figure O The 
integral is performed by splitting into the two regions z > and z < and 
rotating each of the contours to lie on the positive imaginary z axis. This leaves 
six integrals to compute (shown in Figure which combine as follows: 



A+1 y^-\^-\ + ((y2 - A2 - 1)2 - 4A2)i/ 



2A2 7,„i - ((y + A)2- 1)1/2 

/•A-1 // I \\2 _ 1^-1/2 

h + h = +2 / dy — ;^^7, ; J (4.46) 

Jo A2 + 1 - ?;2 _)_ ((y2 ^ \2 _ 1)2 _ 4A2)l/2 

These combine into the simpler integrals 

"2^^* ' -1)"', ("7) 

which are easily evaluated with cosh substitutions. The cutoff b is introduced 
since the separate integrals are divergent. On performing the substitutions we 
find that 

2 <.cosh"^(6+A) 2 /•cosh"^(6+A) 

W{p) ~ lim — / dw coshw— ——r / dw sinh^w 

b^oo A Jcosh-i(A) 2A2 Jcosh-i(6-A) 

_i(A2- 1)1/2 
A 

(L2 _p2)l/2 

= 1 - ^ (4.48) 



so that 



/ d^xda-T{a) = M (l- 
Jp'<P \ 



(j^2_p2)l/2 

L 



(4.49) 

Since the solution is axisymmetric, da-T{a) can only depend on p and z. We 
must therefore have, for p < L, 

da-T{a) = f{p)5{z), (4.50) 

where f{p) is foimd from differentiating (|4.49() : 

M 

= 2.X(L2-p2)i/2 - (4-51) 
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The function f{p) is remarkably simple, given the convoluted route by which 
it is obtained. However, its true significance is not seen until the remaining 
gauge-invariant information has been extracted from Q{a). This information 
resides in the eigenvalues of G{a), the calculation of which introduces further 
complexities. 



4.4 The Einstein tensor 

To calculate the full form of G{a) over the disk, we start with the most gen- 
eral form that G{a) can take consistent with the fact that the Kerr solution is 
axisymmetric. Such a form is defined by, for p < L, 

Q{lo) = S{z){aijo + /3i0 + SiCp + £173) 
g{(j)) = S{z){a2$ + P2'yo + hep + €2^3) 

(4.52) 

Gi^p) = S{z){a3ep + /3370 + Sscj) + £373) 
^(73) = S{z){a4j3 + Pi^i + SiCp + £470), 

where each of the . . . £i are scalar functions of p only. We do not assume that 
g{a) is a symmetric linear function so as to allow for the possibility that the 
matter contains a hidden source of torsion. 

Calculation of each of the terms in g{a) proceeds in the same manner as 
the calculation of the Ricci scalar. The resulting computations are long and 
somewhat tedious, and have been relegated to Appendix ^ The final results 
are that, for p < L, 

g{jo) = 'S{z)j^^^^^^A^{pia + Li,) 
2M 

Gils) = 0. 

These confirm that G{a) is symmetric, so there is no hidden torsion. We see 
immediately that Cp is an eigenvector of G{a), with eigenvalue 2M/{L{L^ — 
p^y^^), and also that there is no momentum flow in the 73 direction, which is 
physically obvious. The structure of the remaining terms is most easily seen by 
introducing the boost factor A via 

, > cosw , . ^ .^ 

tanhA = — (4.54) 

coshw 

and defining the timelike velocity 

V = e^'^'t'jo = coshA7o -I- sinhA^. (4.55) 

(This second use of the symbol v should not be confused with the vector v = IVl 
defined earlier.) With these definitions we see that, over the disk, 

g{v) = 0, (4.56) 
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and 

g[(T^v) = S{z)j^^^^—^^^cT^v. (4.57) 

There is therefore zero energy density in the v direction, with isotropic tension of 
M I {A.-kL{1? — p^Yl"^) in the plane of the disk. This conclusion is gauge invariant, 
since it is based solely on the eigenvalue structure of the Einstein tensor. It is 
truly remarkable that such a simple picture emerges from the complicated set 
of calculations in Appendix 1X1 

The velocity vector v defines the natural rest frame in the region of the disk. 
A second velocity is defined by the timelike Killing vector 

5t = h-i(7o) (4.58) 

(see for details of how Killing vectors are handled within gauge-theory grav- 
ity.) The boost required to move between these velocities therefore provides an 
intrinsic definition of the field velocity in the disk region. In this region the 
function h(a) reduces to the identity, so gt is simply 70. It follows that the 
velocity is given by 

tanliA = cosu = p/L, (4.59) 

and the angular velocity is therefore This is precisely as expected for a 

rigid rotation, which fits in with the observation of |26j that the Kerr solution 
can be viewed as the limiting case of a rigidly-rotating matter distribution. 



4.5 An alternative gauge and the matter ring 

The form of the eigenvectors of Qia) suggest that the more appropriate gauge 
for the study of the Kerr solution is provided by the boost 

R = er^'^fl'^ (4.60) 

so that the new solution is generated by 

W{a) = R{a + Maa-nn)R. (4.61) 

In this gauge the tension lies entirely in the Icr^ eigenplanc, and the character- 
istic bivector of the Riemann tensor becomes 

RcT^R = (4.62) 

which is now a relative spatial vector. In this gauge the 70 frame is the rest 
frame defined by the Weyl and matter tensors, whereas the Killing vectors are 
now swept round. 

The boost H4.6()|l is well-defined everywhere except for the ring singularity 
where the matter is located. This is unproblematic, since the fields are already 
singular there. We know that the integral of the Ricci scalar over the disk gives 

which accounts for the entire contribution to TZ found for integrals outside the 
disk. It follows that the contribution to the stress-energy tensor from the ring 
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singularity must have a vanisliing trace, so that no contribution is made to 
the Ricci scalar. It is also clear that the entire contribution to the angular 
momentum must come from the ring, since the tension is isotropic over the 
disk. From these considerations it is clear that the contribution to G{a) from 
the ring singularity, in the gauge defined by H4.61|l . must be of the form 

4M 

ering(a) = —^(.z) <5(p - L) a- (70 + 0) (70 + (b). (4.64) 

This confirms that the ring current follows a lightlikc trajectory — the natural 
endpoint for collapsing matter with angular momentum. The radius of the orbit, 
L, agrees with the minimum size allowed by special relativity (see exercise 5.6 
of 1231, for example). 



4.6 The physics of the disk 

A natural question is whether the tension field has a simple non-gravitational 
explanation. This is indeed the case. The equations for a special relativistic 
fluid are 

{e + P){v-\/v)/\v ^\/PAv (4.65) 
V-(£w) = -FV-w, (4.66) 

where e is the energy density, P is the pressure and v is the fluid velocity 
(i)^ = 1). For the case of a ring of particles surrounding a rigidly-rotating 
massless membrane under tension we see that (ignoring the factors of S{z)), 
P = P{p) and 

V = coshA7o + sinhA0, tanhA = p/L. (4.67) 

If follows that 

i)-Vt' = sinh^A(^-V0 = - ^ .-, e„, (4.68) 



SO the equation for P is 



This has the solution 



dP p 

^0 



P = 0. (4.69) 



(4.70) 



(L2 -p2)l/2' 

which has precisely the functional form of the tension distribution found above. 
The constant Po is found by requiring that the trace of the stress-energy tensor 
returns M when integrated over the disk. This fixes the tension to M / {AttL{L?' — 
P^)^^^)j precisely as is built into the gravitational fields. Of course, the required 
'light' membrane cannot be made from any known matter. Indeed, the fact that 
the membrane generates a tension while having zero energy density means that 
it violates the weak energy condition. Nevertheless, it is quite remarkable that 
such a simple physical picture holds in a region of such extreme fields. 

In writing the tension as M/ {4t:L{L'^ — p^Y^"^) we are expressing it in terms 
of the radial coordinate p. This coordinate is given physical significance by 
the fact that the h function is the identity over the disk region, so p is the 
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proper distance from the centre of the disk. This gives a simple gauge-invariant 
definition of L as the physical radius of the disk. Furthermore, the form of the 
Riemann tensor (|4.4|l shows that its complex eigenvalues are driven by M/lo^ , 
which is gauge-invariant, and hence a physically- measurable quantity. This 
affords gauge- invariant significance to M and w, and hence to the coordinates 
p and z. It follows that both L and M have simple gauge-invariant definitions, 
without needing to appeal to the asymptotic properties of the solution. 

5 Conclusions 

Many of the significant solutions to the Einstein equations can be represented 
in Kerr-Schild form and gauge-theoretic approach of [2] is well suited to their 
analysis. For all solutions of Kerr-Schild type where the null vector / satisfies 
l-Vl = the Einstein tensor is a total divergence in flat spacetime. The 
structure of the sources generating the fields can therefore be elucidated by 
employing Gauss theorem to transform volume integrals to surface integrals. 
This approach is fully justified within the gauge-theory formulation, since one 
only ever deals with fields defined over a flat spacetime. 

For the case of the Schwarzschild, Reissner-Nordstrom and Vaidya solutions 
the gravitational fields are seen to result from a (5-function point source of mass 
at the origin. For the Reissner-Nordstrom solution the (5-function point source 
is surrounded by a Coulomb field. An unexpected bonus of this approach is 
that the infinite self-energy of the Coulomb field is removed by the gravitational 
field. Similar techniques can be applied to Kinnersley's and Bonnor's work on 
accelerating and radiating masses ISHlEHli as will be discussed elsewhere. 

Applied to more general stationary, vacuum solutions we find that the com- 
plex structure at the heart of vacuum Kerr-Schild fields is the same as the 
natural complex structure inherent in the Weyl tensor through its self-duality 
symmetry. Further algebraic insights are obtained through the use of null vec- 
tors as idempotent elements, simplifying many of the derivations of the vac- 
uum equations. Both of these insights highlight the algebraic advantages of 
the spacetime algebra approach. A further example of this is seen clearly in 
equation (|3.4U|I . which gives a remarkably simple and compact expression for 
the Riemann tensor. 

The application of Gauss' theorem to the Kerr solution reveals some surpris- 
ing features of the singularity. The ring of matter follows a lightlikc trajectory 
and surrounds a disk of tension. The tension distribution over the disk is pre- 
cisely that predicted by special relativity. The correct tension distribution was 
computed by Hamity |14| . though he did not comment on its origin in terms of 
classical relativistic physics. We find no evidence of either the negative surface 
energy density or the superluminal speeds claimed by Isreal Both Hamity 
and Isreal asserted that they used the same results for surface layers in gen- 
eral relativity, but neither gave detailed calculations, so the reason for Isreal's 
disagreement with our result is hard to pin down. 

Almost all trajectories in the Kerr geometry finish up on the disk, rather 
than the ring. Quite what happens when a particle encounters the (5-function 
tension over the disk is unclear and can only really be understood using a 
quantum framework to study the effect of the disk on a wavepacket. (A start 
on such an analysis in made in |3U).) Assuming that all geodesies do terminate 
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on the disk then any non-causal features of the Kerr solution are removed 
which is a physically attractive feature of the picture presented here. The fact 
that the tension membrane violates the weak energy condition raises a further 
interesting question — how can it be formed from collapsing baryonic matter? 
Furthermore, if baryonic matter cannot form the membrane, then what is the 
endpoint of the collapse process? Answers to these questions will only emerge 
when realistic collapse scenarios arc formulated, though these are notoriously 
difficult computations to perform. 

The discussion in this paper implicitly rules out considering any extensions 
to the manifold, such as obtained by converting the Schwarzschild solution to 
Kruskal coordinates. We therefore do not consider distinct universes connected 
by Schwarzschild 'throat', with separate future and past singularities |27[ 
or the maximum analytic extension of the Reissner-Nordstrom geometry with 
infinite ladder of possible 'universes' connected by wormholes |21E2]- In such 
scenarios the applications of Gauss' theorem employed in this paper would not 
be valid. While infinite ladders of connected universes remain popular with 
science fiction writers, there is no reason to believe they could ever form phys- 
ically in any collapse process. The descriptions presented here for both the 
Reissner-Nordstrom and Kerr solutions have a much more plausible physical 
feel to them, even if the final description of the singular region must ultimately 
involve quantum gravity. One final speculation concerns the nature of the mem- 
brane supporting the Kerr ring singularity. This bears a remarkable similarity 
to some of the structures encountered in string theory, and it would be of great 
interest to see if string theory can provide a quantum description of such a 
source. 

A Surface integrals of the Einstein tensor 

From the form of H4.52|l there are 16 scalar functions to find using extensions of 
the technique described in Section ^31 These are evaluated below. 

A.i g{i^) 

For this term we have 

C/(7o) =(5(z)(ai7o + /3i0 + (5iep-hei73), (A.I) 
and ai and ei are computed directly from 

/ ^3x^(70) =2^ /"Vp'(ai(p')7o + ei(p')73)- (A.2) 
Jp'<p Jo 

On converting to a surface integral we obtain 

/'27T /'OO 

d'xg{-io) = P dcj) dzeP-{n{-fo)-jo^{da-n{a))) 

p' <p -^0 J ~oo 

2tt poo / 

2 I o2 



Mp I d(t) dz [ -dpa + {a^ + /3^)ep •71)70 + -90/373 

J-oo V P 



2.AV/ dz(R(7^) + ^i^)7o, (A.3) 

L'^ cosh I 
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which shows that ei = 0. (Here 7"^ refers to the complex field 7 = l/w, and not 
to a reciprocal frame vector. When a 7^-vector is intended the 7 will always 
appear with a subscript.) The final term on the right-hand side of (|A.3|I is the 
integral already performed for the Ricci scalar. For the remaining term we need 



dz 



_^ [p^ + {z + iLYfn- 



where it is again crucial that the correct branch cuts are employed in the eval- 
uation of the contour integral. Substituting this result into ljA.3|) wc find that 

2. J^dp' p'a.ip') ^ AnM (2 - j^r^^ ^ f^j , (A.5) 

and differentiating recovers 

"^ = -L(i2_ ^2)3/2- (A-6) 
For the remaining terms wc first observe that 

d'^xpeP-Gijo) - 27r f dp' p'^ Slip') 
p'<p Jo 

= I d^T(pe''AV)-(r!(7o)-7oA(9a-f^(a))) 

J p'<p 

= 0, (A.7) 

so we must have 5i = 0. To find [3i we need to apply the divergence theorem 
twice: 

27r ['dp' p'^di (pO = / d^x {-p'4>) -5(70) 
Jo Jp'<p 

2ix 



Mp^ 



/ d0 / dz{Ia^)-{-IV(i) 

d^x (/(T3) • (l7(7o) - 70 A (da • f^(a))) 

p'<p 



= 27rMp^(/3(0_) -/3(0+)) + 2M / d^a; (/(T3) • ( V • (a7i)n) 

Jp'<p 

= - (j^2_p2)i/2 + 87rMLW{p), (A.8) 
with W{p) as given by equation H4.48|l . This time, differentiating yields 

Reassuringly, this term vanishes on the axis, as it must do for a valid axially- 
symmctric solution. 
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A.2 ^(73) 

For ^(73) wc can write 

^^(73) =(5(z)(a473 + /?40 + '54ep + e47o). (A.IO) 
This time we find that 

f d^xGijs) ^2tt f dp' p{a4{p')'-f3 + e4{p)jo) 
Jp'<p Jo 

/oo 
dz dp{acrj,-n) 
'OO 

+ Mp-f3 d<P dz{-IcT^)-{^IV{f3<T3-n)) 

Jo J -ex, 

/oo 
dz a sinu 
-00 

= 0. (A.ll) 

The final term vanishes because asinw is an odd function of z. It follows that 
a4 = £4 = 0. The same argument as at equation l)A.7() shows that 64 — 0, and 
for Pi we construct 



27r / V//54(p')= / d'x{-p'4>)-g{j3) 
lp'<p 



= AV/ d4> dz{I(T3)-{~IV{(3a3-n)) 

Jo J -00 

-2 / d^x{I(T3)inh3)-l3/\idaMa))) 
Jp'<p 

/>27r poo 

= ~2Mp d(f> dz {Icr3)-{—a{I(y^)-nn) 

Jo J -00 



AirMp / dz 



°° a sinv cosv 



^ coshw 

= 0, (A.12) 

where we have again used the fact that sin v is an odd function of z. We therefore 
have 5(73) = 0, which is physically reasonable. 

A.3 Gie,), g{^) 
For these terms we write 

g{$) = S{z){a2^ + P2I0 + S2ep + £273) (A.13) 
Giep) 5{z){a3ep + (33lo + ^3^ + £373)- (A.14) 

The calculations are now complicated by the fact that the vector arguments 
are functions of position. To get round this problem we must find equivalent 
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integrals in terms of the fixed 71 and 72 vectors. We first form 

2^ fdp' p'\-2{p') = I dPx{-p'-i^)-g{^) 

= Mp^ d4> dz{~I(T^)-{n{$)-$A{da-n{a))) 

Jo J -00 

+ [ d^x{i(Ti)-{ni-/i)~jiA{da-n{a-))) 
J p'<p 

+ / d''x{Ir72)-{n{-i2)'l2h{da-^{a))) 
J p'<p 

f°° , /3sinhwcosw f°° , asinucosw 

= -2nMp^ dz - h27rAfp / dz ■ 



coshM 7 -00 coshu 

= (A.f5) 

which shows that £2 = 0. A similar calculation confirms that £3 = 0. 

If Q{a) is symmetric then we expect to find P2 = — This is confirmed by 

27T fdp' p'^P2{p')= [ d^Xp'jo-G{$) 

Jo Jp'<p 

= pM dzep-{n{$)~$A{da-n{a))) 
Jo J -00 

d3a;cri.(r!(72)-72A(9„-17(a))) 

p'<p 

d''x(T2-{n{-fi)^jiAida-n{a))) 

p'<p 



f°° acosv 

271 Mp"^ dp dz 2tt Mp dz 

J-00 coshii J_^ 

AnMLp^dp(W{p)/p) - A7:MLW{p) 



_^ coshu 



A similar, though slightly more involved calculation confirms that /Ja = 0. 

For the remaining four functions wc need to consider various combinations 
of integrals. For example, 

7T Tdp' p'{a2{p') + a3{p'))= I d'x7'-e(7i) 
Jo Jp'<p 

2tt 



Mp d(f) dz sm(j){-IcT3)-{-IV{l3(Ti-n)) 
Ja 

np^M 



L 

2ttMp^ 



(/3(0+) - /3(0_)) 

(A.17) 



L(L2_p2)l/2' 



from which we obtain 



MI 2AV /a ion 
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Similarly, 

TT fdp'p'iS.ip') - S,{p')) = f d^Xj'-g{j2) 

= Mpl d(f> sine/) (-/era)- (-/V(/3cr2-n)) 

Jo J ~oo 

= (A.19) 

(because sinhu = over the disk). This result confirms that g{a) is symmetric 
over the disk and therefore that there are no hidden sources of torsion. 
The functions 62 and 63 are obtained from 



J fdp' p'\d2{p') + Ssip')) ^ f d3^/sin,^cos(/)7i-g(7i 
4 Jo Jp'<p 



(A.20) 



which evaluates to 

/>27r poo 

Mp^ d4) I dz sm^4)COfi(t){I(T-i)-{-IV{l3(Ti-n)) 

Jo J -00 

d'x p' COS0 (/era) • (0(71) - 71 A {da-n{a))) 

p'<p 

= -Mp^ i d(p dz siiKpcoscj) a{I An)^ = 0, (A.21) 

Jo J -00 

which shows that 82 = S3 — 0. It follows that Cp is an eigenvector of the 
stress-energy tensor. A similar trick is used to evaluate the final term: 

J /"Vp'^(a2(p')-a3(p'))- / d3x/sin0cos07i-g(72). (A.22) 
4 Jo Jp'<p 

This evaluates to 

^Mp^ d(f) dz sm'^(f)cos(f){I(T3)-{~I'V{f3(T2-n)) 

Jo J -00 

d^x p' cosq>iI(r3)-{ni-f2) - l2A{db -nib))) 
p'<p 

7rMp4 rp 



2L(i2_p2)l/2 



2'n:Mp^W{p) + ATrM f dp' p'W{p'), (A.23) 
Jo 



where W{p) is as defined in equation H4.48|l . We do not need to evaluate the 
final integral since we are only interested in the derivative of the right-hand side. 
This yields 

"^-"^= L(L2-p2)3/2 (A-24) 
which now gives us all of the terms in g{a). 
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